Spontaneous spin textures in dipolar spinor condensates 
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We have mapped out a detailed phase diagram that shows the ground state structure of a spin- 
1 condensate with magnetic dipole-dipole interactions. We show that the interplay between the 
dipolar and the spin-exchange interactions induces a rich variety of quantum phases that exhibit 
spontaneous magnetic ordering in the form of intricate spin textures. 
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Recent achievement of chromium condensate provides 
a platform for cxploring the effects of magnetic dipolar 
interaction in dilute atomic Bose-Einstein condensates 
(BECs) [J. On the theoretical side, early studies on dipo- 
lar condensates 0, concentrated on the scalar atomic 
BECs where the dipole moments are assumed to be po- 
larized by, e. g., external magnetic fields. To unlock the 
magnetic dipole moments, one has to resort to optically 
confmed spinor condensates 0, 0- A widely used and 
very powerful theoretical tool in treating spinor conden- 
sates is the single mode approximation (SMA) in which 
different spin states are assumed to share the same spa- 
tial wave function 0, [?| . An important implication, as 
well as a limitation, of the SMA is the lack of spatially 
varying spin textures, i. e., the spins are uniformly ori- 
ented under the SMA. The SMA is, however, expected 
to be valid only when the total spin-dependent interac- 
tion (spin-exchange and dipolar) is sufficiently weak com- 
pared to the spin-independent interaction j§|, and hence 
cannot cover the whole spectrum of interesting quantum 
spin phenomena. The SMA becomes particularly ques- 
tionable in the presence of dipolar interactions — from 
the studies of ferromagnetic [jj and ferroelectric 01 ma_ 
tcrials, it is known that even a relatively weak dipolar 
interaction may spontaneously induce spatially varying 
dipole moments, a fact which can be attributed to the 
long-range anisotropic nature of the dipolar interaction. 

In the present work, we investigate the ground state 
wave function and spin structure of a dipolar spinor 
condensate by directly minimizing the mean-field energy 
functional without the assumption of the SMA. The main 
result is summarized in Fig. ^ where a phase diagram 
in the parameter space of the relative dipolar strength 
and the trap aspect ratio is present. It can be seen that 
the competition between the long-range dipolar and the 
short-range exchange interactions gives rise to a rich va- 
riety of spin textures. These spin textures emerge spon- 
taneously in the ground state, in contrast to the case of 
non-dipolar spinor condensates, where spin textures can 
only be induced dynamically [llt Il2l Il3| . 

We consider an F = 1 condensate with three magnetic 
sublevels m p = 1,0,-1. In the mean-field treatment, 
the grand-canonical energy functional of the system can 
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where Q (r) = -\/ n a( r ) e i s the macroscopic wave 

function of the mp = a atoms with n a being the density 
and Q being the phase of the component a; /i is the 
chemical potential, n — ^ a n a the total density, V — 
mojQ(x 2 + y 2 + X 2 z 2 )/2 is trapping potential with A being 
the trap aspect ratio, S = J2 a a 4>lFaf34>f3 is the density of 
spin with F being the angular momentum operator, and 
R = r — r'. The collisional interactions include a spin- 
independent part Co = 4irh 2 (a,Q + 2a2)/(3m) and a spin- 
exchange part C2 = / ±irh 2 (a2 — a Q )/(3m) with a f (/ = 0, 2) 
being the s- wave scattering length for two spin-1 atoms 
in the combined symmetric channel of total spin / 
For the two experimentally realized spin-1 condensates 
87 Rb 0| and 23 Na 0], the spin-exchange interactions 
are ferromagnetic (c2 < 0) and anti-ferromagnetic (02 > 
0), respectively. The strength of the dipolar interaction 
is characterized by Cd = ^^gp/^ir) with hq being the 
vacuum magnetic permeability, fiB the Bohr magneton, 
and gp the Lande g-factor. 

The condensate wave function is found by numeri- 
cally minimizing the energy functional E[(j) a , 0*] sub- 
ject to the constraint of fbced total number of atoms 
N using the steepest descent method. In our calcula- 
tion, we fix the value of Co to be positive and choosc 
C2/C0 = —0.01 and 0.03 for the case of the ferromag- 
netic and anti-ferromagnetic spin-exchange interaction, 
respectively. These two ratios correspond to the scatter- 
ing parameters of 87 Rb and 23 Na, respectively 0|. To 
focus on the dipolar effects, we allow Cd to vary and intro- 
duce q = Cdj\ci\ > to measure the relative strength of 
the dipolar interaction. The trap aspect ratio A (which is 
varied between 0.1 and 10 in our calculations) is anothcr 
control knob that allows us to tune the shape of the con- 
densate and hence the effective dipolar interaction 0, [|| . 
Finally, in all numerical results presented in this paper, 
we use N = 10 5 and uj = 2tt x 100Hz. 
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FIG. 1: Phase diagram of dipolar spin-1 condensates of fer- 
romagnetic (a) and anti-ferromagnetic (b) spin-exchange in- 
teraction. s = S/n is the normalized local spin vector. 



Figure ^ summarizes the main result of this paper in 
the form the phase diagram of the dipolar spin-1 conden- 
sate in the q-X parameter space. For the SMA-I, II, and 
III phases, the SMA is found to be valid and the spins 
are uniformly orientcd towards a direction determined by 
the dipolar interaction. In the phases labelled by C, S, 
and P, the spatial wave functions are spin-dependent and 
the system exhibits rich spin textures. 

SMA ground states — The SMA implies that the nor- 
malized local spin vector s = S/n is a constant. For a 
ferromagnetic system in the absence of the dipolar inter- 
action, the vector s has a length of unity and the ground 
state is degenerate about the orientation of s due to the 
SO(3) symmetry possessed by the order parameter 
When the dipolar interaction is present, the SO(3) sym- 
metry is broken and, as shown in Fig. [TJa), the SMA is 
valid only when q is below a A-dependent critical value. 
In the SMA regions, the direction of s is determined by 
the trapping geometry: s is perpendicular to the 2-axis 
for a pancake-shaped (A > 1) condensate (SMA-I) and 
parallel to the z-axis for a cigar-shaped (A < 1) conden- 
sate (SMA-II), reminiscent of a quantum ferromagnet 
with easy-plane and easy-axis anisotropy, respectively. 
This result agrees with the quantum mechanical calcula- 
tion under the SMA @|. 

As shown in Fig. ^b), for an anti-ferromagnetic sys- 
tem, the SMA is again valid for small q. When q is 
below a critical value that is rather insensitive to A, the 
system exhibits a vanishing spin with s = (SMA-III). 
For q above this critical value and a cigar-shaped trap- 
ping potential with A < 0.4, the dipolar interaction may 
overwhelm the spin-exchange interaction and make the 
condensate effectively ferromagnetic, and the system en- 
ters the SMA-II phase with spin oriented along the z-axis. 
We do not find in this case the SMA-I phase with s _L z. 
It is also worth mentioning that transitions between dif- 
ferent SMA phases are all first order. 

Figure ^ shows that the critical dipolar strength at 
which the SMA becomes invalid decreases as A. For a 
ferromagnetic system, q cl - = 0.1 at A = 10, which can 
be reached in 87 Rb |g. The dipolar effects thus become 
more prominent as the condensate becomes more two- 
dimensional (2D) pancake-like. This is consistent with 
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FIG. 2: (Color online). Ground state of the P phase for 
ferromagnetic coupling with A = 2 and q = 1. (a) Density 
of each spin component along the x-axis. The total density 
is denoted by the dashed line. (b) and (c) are the respective 
phase images for cj>i and cf>-i in the z = plane. The phase 
of <j>o is a constant and is not shown here. The phases in the 
z / plane are similar. (d) The vector plot of S in the xy- 
plane with the color map corresponding to the total density. 
Here, as well as in Figs.|3]and2] we adopt a^ = (h/mujo) 1 
and AT/a^o as the units for length and density, respectively. 



the study in solid magnetic materials, where the dipolar 
interaction, normally weak enough to be ignored in bulk 
materials, plays an essential role in stabilizing long-range 
mag netic order in 2D systems, e. g., in magnetic thin films 
[l7f , Our work shows that it is feasible to observe non- 
SMA ground state in a pancake 87 Rb condensate. 

Non-SMA ground states — The SMA regime is anal- 
ogous to the single-domain regime in micromagnetics, 
where the exchange energy dominates over all other spin- 
dependent energy terms and a uniform magnetization 
forms as a result @. As q increases, the SMA eventu- 
ally breaks down. At large q, the dipolar interaction will 
dominate the spin-dependent interaction and the differ- 
ence between the ferromagnetic and anti-ferromagnetic 
spin-exchange interaction becomes insignificant, as can 
be seen from Fig. ^ The non-SMA region is further di- 
vided into three distinct phases P, C and S, which stands 
for pancake, cigar and spherical, respectively (see below). 

Figure |21 shows a typical wave function in the P phase 
which occurs in pancake geometry. The density profiles of 
all spin components are axially symmetric with ni(r) = 
?i_i(r) and hence s z = 0, i. e., the spins are planar and 
lie in the xy-plane. The phases Q take the form 



Q a = w a tp + ip ai 



(1) 



where tp is the azimuthal angle, cp a is a constant phase 
that satisfies 



(f! + (f-! - 2(^o = 0, 



(2) 
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FIG. 3: (Color online). (a) The streamline plot of the spins for 
ferromagnetic coupling in the C phase with A = 0.4 and q = 
1.6. The color map represents |Sj_ | on these cross-sections. 
(b) The vector plot of sii in the za-plane. Same parameters 
as in (a), (c) The same as (b) except for q = 6. 



and w a is the phase winding number with the values 

(wi,w ,W-i) = (-1,0, 1) . 

Therefore the two spin components, <p\ and <fi—i, are in 
vortex states with opposite winding numbers. Due to 
the presence of the non-rotating (j)Q component, the total 
density does not vanish at the vortex core. Such a ground 
state represents a coreless skrymion. It is easy to show 
that the spin vector can be expressed as 



S = 2y/2noni (sm(p, — cos ip, 0). 

As illustrated in Fig. EJd), the spin continuously curls 
around the z-axis to form a vortex such that the system 
possesses a persistent spin current but no net density 
current. The spin vortex state is analogous to the flux- 
closure magnetic states in micromagnetics In the 
vicinity of the z-axis, to reduce the exchange energy, the 
magnitude of the spin gradually decreases and vanishes 
on the z-axis. This differs from the magnetic vortex ob- 
served in nanoseale ferromagnets |18( where, due to the 
conservation of local spin moment, the magnetization in 
the vortex core turns toward the z-axis || . 

In the C phase which oecurs in a eigar-shaped trap, n a 
remains to be cylindrically symmetric and Q can also 
be expressed in the form of Eq. QJ where the phases <p a 
are now z-dependent but still satisfy Eq. (J5J, and the 
winding numbers are now given by 

(w 1 ,w ,w- 1 ) = (0,1,2) . 

The spin vector in the C phase takes the form 
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FIG. 4: (Color online). The spin structure in the S phase 
for anti-ferromagnetic coupling. (a) Spin projection in the 
z = 2.29 plane. A = 0.4 and 5 = 1. (b) Spin projection in the 
y = plane. A = 0.8 and q = 1. (c) The spin structure for 
A = 0.8 and 5=1. (d)-(f) The dipolar strength dependence 
of the toroidal moments for A = 1.2, 0.8, and 0.4, respectively. 



where A = v / 2wo(v / "T + JnZ\), and 5(z) = fo(z) — 
tpi{z) is the spin twisting anglc which is a monotonically 
increasing funetion of z with 5(z — 0) = 0. A typical 
spin texture plot in the C phase is shown in Fig. Eja) 
where one can see that the spins twist around the z- 
axis, tracing out a helical pattern. The total twisting 
angle 8{po) — 8(— oo) increases with q and approaches 
7r at large q. The corresponding structures of the planar 
spin Sii = (S x , S z ) are plotted in Fig. |3fb) and (c), which 
resembles the magnetization of a bar magnet. Note that, 
in the absence of the spin-exchange interaetion (02 = 0), 
the C phase does not exist and is replaced by the SMA-II 
phase with spins uniformly oriented along z. 

In between the P and C phases, is the S phase which 
oecurs when the trapping potential is elose to sphcri- 
cal (A ~ 1). A distinet feature of the S phase is that 
n a becomes non-axisymmetric, signalling the broken of 
the cylindrical symmetry of the spatial wave funetions. 
To further quantify the spin texture, we define the spin 
toroidal moment vector as 01 

G = f dr fr x SMI • 



S = [A s'm(ip + S), —A cos(iy9 + S), 1%% — n_i] , 



For all the SM A phases, we have G = 0. For both the P 
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and C phases, G points along the z-axis. The S phase, 
by contrast, features finite G x and G y . The toroidal mo- 
ments as functions of q is plotted in Fig.0[d)-(f) for three 
different values of A. For small q, G = and the system is 
in the SMA regime. As q exceeds a threshold, G ^ and 
the S phase is reached. The continuity of G suggests, in 
this mean-field calculation, the transitions between SMA 
and non-SMA, and also among non-SMA, phases are of 
second order. In Fig. 0[d), we also see that for q > 0.6, 
the system enters the P phase. Fig. E[a)-(c) show ex- 
amples of the spin structure in the S phase. For this 
particular set of parameters, the condensate displays two 
domains: one has S z > and and the other S z < 0. 

Collapse — Condensate collapse will occur when the 
dipolar interaction strength exceeds a critical value c* d . 
For both ferromagnetic and anti-ferromagnetic coupling, 
we have found that c* d « 0.24co and is nearly independent 
of the trap aspect ration A. This critical value agrees with 
that of a highly elongated dipolar scalar condensate with 
dipole moments polarized along the axial direction 0. 
As is well known, the critical dipolar strength for scalar 
condensate is very sensitive to A [3] • In the spinor system 
studied here, the dipoles are free to rearrange themselves 
to minimize the dipolar interaction, which renders the 
insensitivity of c* d with respect to the trapping geometry. 

Dipole induced spin-orbital coupling — In the SMA 
regime, the total orbital angular momentum L vanishcs, 
and the ground state wave functions cj> a can be taken to 
be real. In the non-SMA regime, the spin and orbital de- 
grees of freedom are intimately coupled and <f) a must be 
described by complex functions. Here the spin (S 2 ) and 
orbital (L 2 ) angular momenta are not separately con- 
served, only the total angular momentum (S + L) 2 is 
conserved. It is not difhcult to show that, in the non- 
SMA regime, if the set of wave functions {<fi a } min- 
imizes the energy functional, then the set does 



not. On the other hand, if we make the transforma- 
tion: (<pi, (fto, 4>-i ) — * — <^o> $*)) then the new set 

still minimizes the energy. The above transformation 
amounts to inverting S and L simultaneously and hence 
preserving the relative orientation between them. This 
inter-connection between S and L is well known in super- 
fluid 3 He and the related phenomenon has been termed 
the spontaneously broken spin-orbit symmetry 19, 20]. 

The dipole induced spin-orbital coupling is recently ex- 
plored theoretically in the dynamics of 52 Cr condensate 
and is predicted to manifcst in the Einstein-de Haas ef- 
fect [21|. 52 Cr atom features a spin-3 ground state j2^| 
whose short-range collisional interaction is characterized 
by four scattering lengths and not all of them are accu- 
rately known. Although we have studied, in this work, a 
relatively simpler spin-1 system, we expect much of the 
essential physics can be applied to higher spin systems. 

In conclusion, we have provided a detailed phase dia- 
gram which shows the ground state structure of a dipo- 
lar spin-1 condensate beyond the SMA. In this system, 
the spin and orbital degrees of freedom are intimately 
coupled together. The interplay between the long-range 
dipolar and the short-range exchange interactions gives 
rise to a variety of quantum phases characterized by dis- 
tinctive spin textures. As such, dipolar spinor conden- 
sates represent an intriguing quantum magnetic system 
whose properties are highly tunable. The work here as- 
sumes zero magnetic field, an assumption valid for 87 Rb 
and 52 Cr when the external magnetic field strength is less 
than 0.1 mG @. This constraint poses an experimental 
challenge but is definitely within the reach with current 
technology. 
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